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We investigate transport properties of gate-all-around Si nanowires using non-equilibrium Green's 
function technique. By taking into account of the ionized impurity scattering we calculate Green's 
functions self-consistently and examine the effects of ionized impurity scattering on electron den- 
sities and currents. For nano-scale Si wires, it is found that, due to the impurity scattering, the 
local density of state profiles loose it's interference oscillations as well as is broaden and shifted. In 
addition, the impurity scattering gives rise to a different transconductance as functions of tempera- 
ture and impurity scattering strength when compared with the transconductance without impurity 
scattering. 
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I. INTRODUCTION 

The study of ballistic electron transport in nano- 
devices has been an interesting field of research. 
Recently, a Si nanowire with a length comparable to the 
de Broglie wavelength of carriers is realized by advanced 
nanofabrication technique. [5, 6] The cross-sectional area 
of Si nanowires was designed to show well-separated 
transverse modes and electrons confined to the wire are 
expected to suffer from a minimal amount of impurity 
scattering. These properties make the Si nanowires good 
candidates for the study of ballistic quantum transport. 
In addition, the potential distribution within the wire 
can be controllable by a metallic gate around the wire. 
This provides additional degree-of freedom on currents 
through the device and one would expect that the basic 
transistor action is possible for a Si nanowire. As a re- 
sult, the gate-all-around Si wire may shed the light on 
one-dimensional structures for future transistor applica- 
tions. 

It is desirable experimentally to make the Si wires as 
intrinsic as possible. However, to populate the wires with 
carriers, it is necessary to define source and drain regions 
where ionized dopants are placed. These dopants scatter 
free carriers and the elastic impurity scattering cannot be 
avoided in those regions. Thus, in order to understand 
transport in the wires, a quantitative treatment of the 
ionized impurity scattering will be important. Several 
theoretical works were done to investigate the effects of 
ionized-impurity scattering on one-dimensional electron 
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gas, and revealed their effects on the electronic struc- 
ture. Most of these studies were for uniformly doped 
or remote-impurity systems 0, Q and adopted empirical 
models based on the so-called Biittiker probes for simu- 
lating the device. [1, [13] The empirical methods are ap- 
pealing due to relatively simple implementation but the 
methods often require parameters that need to be ad- 
justed using more rigorous calculations or values from 
experiments. 

In this work, we take into account the ionized impu- 
rity scattering in simulating the gate-all-around nanowire 
using non-equilibrium Green's function approach. By 
averaging the Green's function over impurity configu- 
rations and expanding the arising term perturbatively, 
we treat the impurity scattering within a self-consistent 
Born approximation and apply the formula to the Si 
nanowire as realized in Ref. [6[. Since the impurity- 
scattering strength is a single parameter for the sys- 
tem, the method provides the first-principle approach to 
understand current-voltage characteristics and compare 
them with the experimental results. 



II. CALCULATION METHOD 

A. Hamiltonian 

To see the effects of the impurity scattering clearly, 
we consider a simple geometry of a quantum wire as in 
Fig. [TJ An infinitely-long cylindrical Si wire consists 
of intrinsic channel and heavily doped source and drain 
regions. A metallic gate extended over a length of Lq is 
rolled round the intrinsic region and they are separated 
from each other by a Si02 layer with a width <ox- For 
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FIG. 1: We plot a schematic diagram of a cylindrical Si wire 
simulated in this work which is oriented along the [001] di- 
rection. The Si wire surrounded with the gate is assumed to 
be intrinsic and separate the source and drain regions where 
ionized dopants are distributed. 



simplicity, we assume that the Si wire is grown along the 
crystal [001]-axis(chosen as the z direction in the figure) 
and the doping profile of Njj (r) in the source and drain 
regions is symmetric about the 2;-axis so that we can 
utilize the circular symmetry. 

Then, electrons in the Si wire are governed by the 
effective-mass Hamiltonian which is given by 
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Above Hamiltonian describes electrons in six different 
valleys depending on their effective masses. For instance, 
if nix — — 0.19m, transverse mass, and rriy = 0.95m, 
longitudinal mass of Si, the Hamiltonian represents elec- 
trons in the [010]-valley, etc. Here, U{r) is the macro- 
scopic potential energy resulted from both band discon- 
tinuity among the materials, and the Coulomb contribu- 
tion from external charges. The Coulomb part is deter- 
mined by the Poisson's equation. 



^y^U{r) = —{NDir)^nei{r)} (2) 

when we know the electron distribution 7i£/(r). U^mpi^) 
describes the impurity potential energy from the ionized 
dopants. In this work, we assume that the impurity 
potentials are short-ranged but still vary slowly in the 
atomic scale. As a result, different valley modes are not 
coupled by the impurity potential and can be solved in- 
dependently. 

Since the device has the circular symmetry, it is con- 
venient to express the Hamiltonian in terms of the basis 
diagonalizing the radial motion. We choose the basis sat- 
isfying the following Schrodinger equation. 



1 d^ 



1 d^ 



^nix dx^ niy dy"^ 



-Ub{p) 



\xi)^^i\xi) (3) 



where p is radial coordinates {x,y) and Ul{p) = 
U(p,zLoo) is a potential energy aX z — ±oo, i.e., in the 
deep source and drain regions. Then, we expand the field 
operator V;(r) as. 



(4) 
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where we discretize the longitudinal coordinates with a 
spacing of a and tpizm) is tight-binding basis at the m— th 
node {zm = mo, m = — oo, ..oo). 

Using Eq. ([4]) and a finite difference approximation, 
one can express the Hamiltonian of Eq. ([1]) as, 
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(5) 



Here, the first term describes motion along the longitudi- 
nal direction for each transverse mode and it's elements 
are given by, 



H/m:;'m' = [5m,m'{^lm + 2t/J-) — i// ((5m,m' + l + <5i) 



with e„ii^ei + {xi\U{p,Zm)-UL{p)\xi) and the hopping 
energy of tu — jlmzo} (hereafter, we use bold charac- 
ters to denote a matrix displayed on the basis {x/V'm})- 
The V matrix in Eq. accounts for the deviated po- 
tential distribution from that of deep source and drain 
regions. As a result, it gives rise to the hybridization 
among transverse modes as, 

V;m:/'m' = Wl\U {p, Z^)-U l{p) \Xl){^ ^ ^l.v)^mnA^) 

The last term in Eq. ([5]) is a contribution from the im- 
purity potential. 



B. Impurity-averaged Green's function 

Now we formulate non-equilibrium Green's functions 
for the Hamiltonian of Eq. . In order to take into ac- 
count the impurity scattering, we consider a number of 
impurity configurations rather than a particular distri- 
bution, and average the Green's functions over the con- 
figurations. For this we adopt the Schwinger-Keldysh 
technique, [m According to the scheme, the impurity av- 
erage gives rise to the quadratic interaction in the ac- 
tion, and we expand it perturbatively to obtain the one- 
particle irreducible self-energy S™^'. Here, we restrict 
our attention to the first order diagram and treat it self- 
consistently, which is referred to as the self-consistent 
Born approximation. [X, ^, JJi, [l3| 

The impurity-averaged Green's function G can be ob- 
tained through the Dyson's equation. 



(7) 



where g is the impurity-free Green's function (in fact, the 
bold characters in this case represent enlarged matrices 
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for taking into account the Keldysh space. However, we 
keep the notation in the meanwhile because it recovers an 
original size when we specify it's components explicitly in 
the Keldysh space). The corresponding self-energy from 
the impurity scattering depends on it's Green's functions 
again through the relation, 



with 
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2\(XiV'm|?7imp(r) I XlAnii) 
(XisV'ma |t^«mp(r')IXi''0m') / -(9) 
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Here, (• • • )av denotes a configuration average. We model 
fluctuating impurity potentials with a (5-correlated func- 
tion considering the short-ranged form; 

{U^.mpir)U^„,p{r'))av = nD{r)ul S{r - r'). (10) 

Here, n£)(r) — ND{r)/No is a normalized doping pro- 
file with respect to the atomic density A^o of Si. And 
the impurity potential strength is expressed with the im- 
purity potential amplitude of uq and a screening length 
Is = 4A, which is approximately equal to the Tomas- 
Fermi screening length in the bulk Si at carrier density 
of 1 X 10^°/cm'^. Accordingly, the expansion coefhcient 
in Eq. ^ becomes, 
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It is noted that the short-ranged potential is diagonal 
for longitudinal basis {ipm} but not for transverse modes 
{xi}- This means that transverse modes are mixed to 
each other through the impurity scattering. 

For a given in order to solve the Dyson equa- 

tion of Eq. ([7]) , we should take care of open-boundaries 
in our problem, i.e., the infinite number of nodes along 
the longitudinal direction (m — —00, ...,00). For this, 
we follow the conventional approach where the device 
is partitioned into the system being in non-equilibrium 
and reservoirs. [1^ Since the source and drain regions 
are extended semi-infinitely, we confine our attention to 
the portion of the system near the gate where physical 
properties are thought to be deviated from those of deep 
source and drain regions. We designate the portion by 
longitudinal indices m = (0, 1, M— 1). Thus, nodes for 
TO < (to > M) represent the source (the drain) being 
in equilibrium with the chemical potential /is (pd)- 

In the source and drain reservoirs, we assume that 
the self-energy S""p is independent of longitudinal co- 
ordinates TO because they are sufficiently far from the 
gate region where the potential distribution is uneven. 
Within this assumption, the Schrodinger equation is 
easily solved and equilibrium Green's functions G{E) 



with corresponding self-energies are calculated straight- 
forwardly. In the Appendix, we illustrate their simple 
expressions. 

Now, we focus on the device region, i.e., nodes rang- 
ing < m < Af where one expects a non-equilibrium 
situation for different chemical potentials of fis and fiD- 
The Green's functions are obtained by truncating the 
matrix equation of Eq. ([7]) within longitudinal indices 
of < TO < M. Instead, the truncation introduces an 
additional self-energy S to the Dyson equation owing to 
the coupling of the source and drains, and a total self 
energy becomes S = S + S'™p. Here, the self-energy 
t,{E) reads, 

^lm:l'm'{E) — t'j^Smm' '5m,oG;(_ i) (_ i) (i?) | p=^g 

+^m,M-lGlM:l'M{E)\^=fj,o (12) 

where the subscripts oi ^ — hs.d denote that each equi- 
librium Green's function is determined by different chem- 
ical potentials of ps = Po — ^Vs and = Po — eVo 
accounting for applied voltages, Vg and Vd at each reser- 
voir, respectively. 

Solutions of the Dyson equation are obtained by in- 
verting the matrix equation Eq. ([7]). Firstly, it's retarded 
component is calculated as. 



G«(£;) = [(g«) 



(13) 



Here, g^{E) = [SI — H — v] is the free-particle Green's 
function and S^(£;) = t^{E) + 'S'"'P'^{E) is a retarded 
component of the self-energy. Detailed form of T:^{E) is 
given in the Appendix. Whereas, the term of S*''"P'^(£') 
depends on diagonal components of it's own Green's func- 
tion, as indicated by Eq. ([5]). Thus, we should solve the 
above matrix equation self-consistently. 

With the obtained G''^ and it's Hermitian conjugate 
G"^, the Keldysh components of the Green's function and 
the self-energy become 

G^{E) = G'^iE)'S^{E)G^{E) (14) 



and 



^^{E) ^ t^{E) + i:'"'P'^{E), (15) 



respectively. According to Eq. the self-energy con- 

tributed from the the source and drain coupling is ob- 
tained as, 
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2kBT 
Pd~ 



2kBT 



) (16) 



with t:^{E) = 'E'^{E) ~ S^(£;), the correlated compo- 
nent of the self-energy. However, for the Keldysh compo- 
nent of the impurity- induced self-energy S""P'^(i?) the 
result is not given in a closed form and should be calcu- 
lated self-consistently as in the case of the retarded one 
via Eqs. © and (O). 
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C. Electron density and current 



D. Approximations 



The ensemble average of n;„i = {b\„^bim) gives local 
electron density of the device and, consequently, the elec- 
tron density distribution in Eq. ^ becomes Ueiir) = 
J2 im '^irnXi (p)^^ {zm) ■ ^Fvom the generating functional 
technique as in Rcf. 15], one can express the local elec- 
tron density with the calculated Green's functions. The 
result reads, 



Za 



tr 



dE Gg,.,„,iE) 



dE 



iFD{E)Ti{lm : E) 



(17) 



Here, in the second line we use the functional form of 
Fermi-Dirac distribution fp£){E) and the density-of-state 
D(/to : E) for the resemblance with equilibrium results. 
Since the device is in non-equilibrium condition, two 
functions are given in a matrix form; the Fermi-Dirac 
distribution matrix is defined by. 



(18) 



while, using Eq. ([H]) . the density-of-states matrix at the 
node m and transverse mode I, is expressed by. 



ZTra 



(19) 



Here, gsv = 4 is the spin- valley degeneracy, S*-^ = 
— S'^, and 1/m is a matrix whose elements are non- 
zero only at the Zm-th diagonal position. When the im- 
purity scattering is absent, tpo becomes the well-known 
results as in Ref. [13, [l3| , where non-zero elements are 
only at m = and m = M — 1 nodes and are equal to 
the Fermi-Dirac distribution characterized hy fj,s and/ic 
respectively. However, due to the impurity scattering of 
elements oifpjj are deviated from the Fermi-Dirac 
distribution function in general. 

Currents flowing through the device is defined by time- 
derivatives of total charge at nodes m = —1 or m = M . 
Then, through the Heisenberg equation of motion, one 
can find that the currents becomes, 



Ids 



2Trh 
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trSR 



dE 



/oo 
dE iFDiE)T,n{E) 
-00 



(20) 



where, by to = or i\f — 1, the expression means 
currents at the source or the drain, respectively, and 
Im = X); '^im- In the second line of the above equation, 
we define the transmission matrix by. 



Prior to numerical calculations, let us first look at the 
approximations used. Firstly, we consider a finite number 
N of transverse modes. Then, the solution of Eq. (|13p is 
obtained by inverting a (NM) x (NM) matrix iteratively. 
However, this scheme demands the huge computational 
cost because the matrix size is large and is deviated from 
the tridiagonal form due to off-diagonal elements of the 
self energy S and the Hamiltonian v. 

As an approximation, we consider leading terms in 
Green's functions to emphasize mainly the effects of the 
impurity scattering. This is equivalent to consider the 
diagonal components of the Green's functions for trans- 
verse modes. Namely, the coupling of different transverse 
modes in the self energy S and the Hamiltonian matrix v 
are neglected. As indicated in Ref. [l3|, if the potential 
energy U{r) is a slowly-varying function along the radial 
direction at any node to the Hamiltonian matrix v be- 
comes small and the approximation is well justified. As 
for the self-energy, leading terms in the Green's functions 
are obtained by writing overlap functions of Eq. (jllD as, 

SlmVm-.hmhm Kh-TT'iXl II Xhip) noip^Zm) \ Xl') 

la 

^5i,ijiv^{xi\\ XiAp) \'nD{p,Zm) lx/022) 

and, therefore, the self-energy of Eq. ([5]) becomes diag- 
onal for transverse modes. However, the approximation 
of Eq. ((22)) still couples transverse modes non-trivially 
because each diagonal component of the self-energy de- 
pends on others. 

Another approximation is made in the Keldysh com- 
ponent of the impurity self-energy s*™^''^. After various 
numerical calculations, we find that S""^^^ is well rep- 
resented by; 



''Ini-.l'm 



,{E)=Su,Sr 

i tanh(^; 
1 tanh(|^ 



for TO < M/2 
for TO > M/2 



(23) 



where a node to = M/2 is the middle point in the in- 
trinsic Si wire. This indicates that particles at the nodes 
near the source(drain) have still the chemical potential 
Ps (pd), not an intermediate value between fis a-nd fijj, 
even after suffering from scattering. We attribute this 
result to a particular potential distribution in the device 
of a source-to-channel barrier, which prevents particles 
with different chemical potentials from mixing. 



III. RESULTS AND DISCUSSIONS 



Tm — 3si>5] (^l.,„G S S !,,„ — G j (21) 

In the case of free inipurities, this form also recovers the 
previous results. [13, [3| 



In this section, we numerically illustrate solutions of 
the non-equilibrium Green's functions suffering ionized 
impurity scattering and related transport properties. We 
consider a typical case of the device structure which can 
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be realized experimentally. As shown in Fig. [Tl the 
source and drain regions are doped at 10^° /cm'^ and 
there is no gate-to-source and -drain overlaps to consti- 
tute nearly abrupt junctions with the intrinsic channel. 
The source and drain extensions are 15nm and the gate 
length Lg is 20nm, so that a total device length simulated 
is 50nm. By choosing a node spacing of a = 0.25nm, we 
have the number of 200 nodes along the wire. In order 
to highlight quantum effects, we choose a small radius 
(3nm) of the wire which exhibits 3 mode occupancies at a 
zero temperature. However, to include thermally excited 
particles as well as the mode coupling from the impurity 
scattering, 20 transverse states are incorporated. The 
gate oxide layer has a thickness of 2.5nm and is treated 
as an infinite potential barrier for electrons. Due to this, 
wavefunctions at the interface between the Si wire and 
the oxide are assumed to be zero in all of our simulation. 

The Poisson's equation is solved in the cylindrical 
coordinates with Dirchlet boundary conditions at the 
gate-oxide interface, otherwise, with Neumann condi- 
tions. For a rapid convergence of solutions, we use the 
Newton- Rhapson method for the Gummcl form of exter- 
nal charges, [lg To model a gate material, we choose a 
work function of 4.56ey, approximately for TiN. 

In Fig. [21 we show calculated electronic subbands 
of each level and local particle density along the wire, 
and compare the results with and without the impu- 
rity scattering in (a) and (5), respectively(VG = 0.6F 
and Vd = O.IF). The subband bottoms(dotted lines) 
reflect the calculated self-consistent potentials in which 
electrons at each levels feel at a node m. Regardless of the 
impurity scattering, they exhibit source-channel barriers. 
Since a high gate voltage lowers the energy barriers, the 
basic transistor action is achieved by controlling these 
barriers. [l3| 

The energy-resolved particle density is plotted in a 
gray scale; a darker area in the figure represents higher 
density. In the impurity- free case of (a), since there 
is no momentum relaxation, states injected from the 
drain (source) end of the device undergo reflections and 
interfere strongly to the right(left) of the source-to- 
channel barrier. This interference results in coherent os- 
cillations in the particle density as seen in Fig. [ll-(a). As 
a function of energy, it is found that the local particle 
density far from the source-channel barrier shows sharp 
peaks like 1/\/E at every onset of subbands, reminiscence 
of one-dimensional density of states. 

If one turns on impurity scattering, phase information 
of the electrons within the device is randomized and the 
energy levels are renormalized. Above all, this makes the 
interference oscillations washed out in the local particle 
density as shown in Fig. [l]-(b). In addition, electronic 
states are shifted and broaden, so that the most electrons 
are found below subband bottoms and it's occupation has 
no longer 1 / -s/E-dependence, but a monotonically vary- 
ing function(the abrupt change of darkness along the en- 
ergy direction comes from a different valley state). In 
both cases of the impurity scattering, one can see that 
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FIG. 2: For a cylindrical Si wire with a 20nm gate length, 
we plot local particle density nim along the wire for im- 
purity scattering strengths of = and 39(eV) in (a) and 
(6), respectively, at Vgs = 0.6V, Vd = O.IV, Vs = OV, and 
T — 300K. For qualitative comparison, we display higher den- 
sity with a darker color. Dotted lines describes the effective 
potential energy of each subband before they are renormalized 
by impurity scattering. 



electrons in the source and drain regions are well sep- 
arated by the source-channel barriers from each other. 
Due to this, the approximation of Eq. ((23)) is justified 
with good accuracy. 

In order to examine the electronic transport of the de- 
vice, we calculate channel currents Ids versus a gate volt- 
age Vg at a small source-drain bias, and plot results in 
Fig. [3]- (a) and (b), respectively, with and without impu- 
rity scattering for several temperatures. Under this con- 
dition, currents exhibit rapidly increasing behavior as a 
gate voltage becomes larger. This shows the basic oper- 
ation of a transistor as indicated in the previous section; 
the channel current turns on by lowing the source-channel 
barrier when a gate voltage is higher than a certain value, 
called a threshold voltage Vth- 

By comparing Figs. [5]-(a) and (b) at a given tempera- 
ture, one can find that the presence of impurities reduces 
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FIG. 3: We compare calculated Ids-Vg results at tem- 
peratures of 300K(solid), 200K(dotted), lOOK(daslied), and 
50K(dot-dashed lines), respectively, for impurity scattering 
strengths of «o = in (a) and 39eV in (6). Here, we assume 
a small source-drain bias of 0.021^. 



the currents significantly even though electrons in both 
cases are expected to move ballistically in the intrinsic 
gate region. This indicates that transport through the 
Si wire largely depends on the electronic structure of the 
source and drain regions. 

As inspired by flat subbands in the figures, the poten- 
tial drops across the intrinsic regions are nearly invari- 
ant to the impurity scattering strength. Thus, it is rea- 
sonable to assume that the suppressed currents do not 
come from the Fermi-Dirac matrix of Eq. (|20p which 
crucially depends on the potential drop, but mainly from 
a reduced transmission coefficient of Eq. pT|) . One of 
possible explanations for this is that electrons injected 
from the source are partially refiected from impurities in 
the source extension in addition to that from the source- 
channel barriers and, thus electrons tunnel the source- 
channel barrier at rare intervals. This type of the reduc- 
tion for the transmission coefficient is also encountered in 
problems of tunneling in dissipative environments, [l9j 
According to the theories, when environments of the de- 
vice become more dissipative, carriers are harder to tun- 
nel the barriers because more energies should be trans- 
ferred to the environment. 

As a function of a temperature, curves are shifted with 
wholly similar shapes and slightly different slopes in both 



cases of the impurity scattering. Two points are note- 
worthy. Firstly, the threshold voltage is shifted to a 
higher value as a temperature is lowered. This is easily 
understood because as the temperature decreases, avail- 
able electrons to overcome source-to-channel barrier ther- 
mally are reduced and then more potential energy should 
be supplied electrostatically to turn on currents. Sec- 
ondly, we look at the slopes of the Ids-Vg curves. In con- 
ventional MOSFETs, they are related to a channel mo- 
bility fj,m via a relation of Ids oc UmiVc - Vth){VD - Vs). 
As seen in the figures, our results show linear behavior 
in some range of gate voltages. Therefore, we may un- 
derstand the slopes to be proportional to the mobility 
of electrons in the device. For detailed comparison, we 
define the conductance by 



a(T) 



dlos 



dVa 



(24) 



Vg=0.65V 



known as the transconductance in MOSFETs. 

Calculated conductance is summarized in Fig. [5] as 
functions of impurity scattering strength and tempera- 
ture. In Fig. HI- (a) we compare the conductance with 
increasing impurity scattering strength for two tempera- 
tures. It is noted that the conductance decreases mono- 
tonically when the impurity scattering strength becomes 
larger at both temperatures and, consequently, sup- 
pressed mobilities are expected. In Fig. [ll-(b) we plot 
the temperature dependence of the conductance for var- 
ious impurity scattering strengths. For a bulk material, 
it is well known that the mobility resulted from impurity 
scattering is proportional to T^/^ to the first order(dotted 
line in the figure). [20] In the case of a two-dimensional 
system, the ionized impurity scattering (for instance, in 
a quantum well with a 5-doping) is enhanced due to the 
increased overlap of the ionized impurity with electron 
wavefunctions and the mobility decreases nearly expo- 
nentially when a temperature is lowered(dashed line).[2l| 
In our case of a quasi-one-dimensional system, the con- 
ductance shows different temperature dependences from 
those of higher-dimensional ones; the conductance of the 
Si wire interpolates from linearly increasing behavior of 
the impurity-free case to the exponentially decaying de- 
pendence of a strong impurity scattering as a function of 
scattering strength. Curves shown in Fig. [ll-(b) do not 
provide a definitive comparison of ionized-impurity scat- 
tering among three different dimensional systems because 
each system has different doping profiles and concentra- 
tions. Despite of this, it is interesting to note that the 
ionized impurity scattering becomes less temperature- 
dependence when the system has a lower dimension. 



IV. SUMMARY 

In summary, we study transport through a gate-all- 
around Si wire in the ballistic regime by considering 
the ionized impurity scattering. Using the Schwinger- 
Keldysh approach, we include the impurity scatter- 
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ing within the self-consistent Born approximation and 
present expressions for electron density and currents 
in terms of non-equilibrium Green's functions and self- 
energies. By simulating a typical case of a Si wire, we 
compare electron densities and channel currents for zero- 
and strong-impurity scattering strengths. In the case of 
the strong impurity scattering, we find that the local par- 
ticle density profiles are shifted and broaden to result 
in suppressed currents compared to the zero-impurity 
scattering case, and the oscillating interference pattern 
vanishes. Calculated currents and conductances are also 
presented as functions of temperature and the impurity 
scattering strength. It is found that the conductance of a 
Si wire exhibits various behavior by decreasing tempera- 
ture, which interpolate from a linear increasing function 
at a zero scattering to an exponentially decreasing func- 
tion for the strong scattering case. However, in this work, 
we do not include other inelastic scattering process such 
as acoustic and optical phonon scattering which will be 
occurred in real devices. Therefore, our results show the 
effects of the ionized impurity scattering alone on elec- 
tronic transport through a Si wire. 
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FIG. 4: In (a), we plot calculated conductances(symbols) as 
a function of impurity scattering potential at two different 
temperatures of 300K and lOOK, rcspcctivcly(VD = 0.02^). 
In (6), calculated conductances arc jjlottcd as a function 
of temperature for given impurity scattering potential of 
Mo = O(circles), 23ey (crosses), and 39el^ (triangles), respec- 
tively. Solid lines are just guide to the eye. To emphasize 
their temperature dependence we normalize them with values 
at 300K and superimpose the lines of gi ^^Cr/sooK-i) (-(jashed) 

and (T/300K)3/2 (dotted). 
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with 

' AtH 

and, according to Eqs. ^ and (P^ . the self-energy is 
proportional to the diagonal component of the Green's 
functions like 



APPENDIX A: GREEN'S FUNCTIONS OF A 
UNIFORMLY-DOPED WIRE 



'Sw 



h 



Jlrnlrn-.lirn-.lim^l-^mdim 



(A3) 



In this Appendix, we illustrate Green's functions for an 
infinitely long Si wire which is doped uniformly. Eigen- 
states are plain waves whose wavelength is determined by 
periodic boundary conditions. Since the wire is transla- 
tional invariant, the self-energy in Eq. ([5]) is independent 
of a longitudinal position. Then, the retarded component 
of the Green's function can be derived as, 



Ma/ 2 



lm:l'm'\ / 71 //- / , ^ 
Mo — , hj 

k=-Mo/2 



^27rik{m — 'm') /Mq 



^Ik 



•^irnp.R , 



(Al) 



where Mq is the number of a longitudinal node and 
f-ik = e/ + 2tjy[l — cos{2iTk / Mq)] with an eigenenergy 
of the Z— th transverse mode. In the limit of a large Mq, 
diagonal components of the Green's function reads; 



1 sign(y;(^)-l/2) 
4tH ^yfiE)-yiiE) 



(A2) 



Thus, the Green's functions are obtained by solving 
Eqs. (jA2p and (|A3|) self-consistently. On the other 
hand, the chemical potential fiQ of the uniformly doped 
wire can be found from the particle density of n/,„ = 
— gst,3G^j.j„j(£^)/7ra together with the Poisson's equa- 
tion. 

The self-energy of Eq. caused by the coupling of 
the device to the source and drain regions is obtained by 
solving the uniformly doped wire with vanishing bound- 
ary conditions. In the similar way to Eq. (|A2p . it is given 

by, 



Su'tH[2yi{E)-l] 



1- Wi 



(A4) 



